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Abstract. We discuss the quantum dynamics of the PU oscillator, i.e. the system with the 
Lagrangian 



L = - [q 2 - (fij + nl)q 2 + 0^2?] (+ nonlinear terms). 



(1) 



When ill 7^ O2, the free PU oscillator has a pure point spectrum that is dense everywhere. 
When fix = r?2, the spectrum is continuous, E g {—00,00}. The spectrum is not bounded 
from below, but that is not disastrous as the Hamiltonian is Hermitian and the evolution 
operator is unitary. Gcnerically, the inclusion of interaction terms breaks unitarity, but 
in some special cases unitarity is preserved. We discuss also the nonstandard realization 
of the PU oscillator suggested by Bender and Mannheim, where the spectrum of the free 
Hamiltonian is positive definite, but wave functions grow exponentially for large real values 
of canonical coordinates. The free nonstandard PU oscillator is unitary at f2i 7^ n 2 , but 
unitarity is broken in the equal frequencies limit. 
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1 Motivation 

Mechanical and field-theory systems involving higher derivatives in the Lagrangian attracted 
(or, better to say, reattracted) recently some attention. In [TJ [2] we put forward arguments that 
the Theory of Everything may represent a variant of supersymmetric higher-derivative theory 
living in flat higher-dimensional space. Our Universe is associated then with a 3-brane classical 
solution in this theory (a kind of soap bubble embedded in the flat higher-dimensional bulk), 
while gravity has the status of effective theory in the brane world-volume. 

The first paper where such theories were considered dates back to 1950 [3]. It was shown 
there that higher-derivative theories are in many cases ghost-ridden, which may break unitarity 
and render theory meaningless. Recently, it was understood, however, that there are higher- 
derivative systems where the ghosts are "benign" such that unitarity is preserved. An example 
of nontrivial supersymmetric higher-derivative system of such benign variety was constructed 
and studied in [4j. 

In this note, we concentrate on the dynamics of the system ([1]) (it is the simplest example 
of higher-derivative mechanical systems introduced and studied first in the paper [3], which is 
known by the name Pais- Uhlenbeck oscillator). Recently, the interest to this problem has been 



*This paper is a contribution to the Proceedings of the Vllth Workshop "Quantum Physics with Non- 
Hermitian Operators" (June 29 - July 11, 2008, Benasque, Spain). The full collection is available at 
|http: //www.emis.de /journals/SIGMA/PHHQP2008.htmT1 

i On leave of absence from ITEP, Moscow, Russia. 
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revived [21 [5j [9] . In spite of that many salient features of the PU oscillator dynamics were 
revealed in previous studies, a certain confusion persists now in the literature. This especially 
concerns the interpretation of the results obtained. Thus, we found it useful to write a kind of 
mini-review including all relevant previous results with inaccuracies corrected, as well as some 
new remarks. There are basically three new observations: 

• The statement of [5] that the Hamiltonian of the oscillator in the equal frequency limit 
represents a set of Jordan blocks of infinite dimension is compatible with the statement 
that the spectrum of such system is continuous if the limit £l\ — > O2 is defined in a natural 
way [31 [6] . We note that the wave functions of continuous spectrum can be represented 
as superpositions of non-stationary solutions to the Schrodinger equation (non-stationary 
solutions are characteristic for Jordan block Hamiltonians.) 

• In contrast to what we suggested earlier [5], the spectrum stays unbounded from below 
also when nonlinear terms in equation ([I]) are included. In most cases, the latter bring 
about the quantum collapse phenomenon and breaking of unitarity. Unitarity survives, 
however, when interaction has a certain special form. 

• In contrast to what was written in [9], the free equal- frequency PU oscillator in the non- 
standard realization is not unitary in the usual sense of this word. 



2 Hamiltonian and its spectrum 

The canonical Hamiltonian corresponding to the Lagrangian ([1]) can be obtained using Ostro- 
gradsky's method [ID] . Due to the presence of extra derivatives, the phase space involves besides 
(q,p q ) an extra canonical coordinate 

x = q 

and the corresponding canonical momentum p x . The canonical Hamiltonian is 

H = p q x + q+ V 2 2 > (2) 

Indeed, if writing the canonical Hamilton equations of motion and excluding p q , x, p x , we arrive 
at the equation 

g (4) + (nj + nj)q + n\nlq = o, 

which can on the other hand be directly derived from the Lagrangian (jTJ) . 

Consider first the case of inequal frequencies and assume for definiteness f^i > ^2- The 
Hamiltonian ([2D can then be brought into diagonal form by the following canonical transforma- 
tion 

_ 1 OiX 2 - Pi _ ShXi - P 2 

Q ~ ^ly/Q-i-n^ X ~ - ' 

\J — \J ^1 — ^2 

The inverse of it is 

l p q + Q 2 x Px + Mjq 
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We obtain 



The spectrum of this Hamiltonian is 

^nm=(n + 7)^i-(m+o) fi 2, n,m = 0,1,2, ... . (6) 



2 / V 2 
The eigenfunctions of the original Hamiltonian ([2]) are [5] 



where A = Oi — O2 and 



m < n, 



_ ^ / zA \ m!(n-m)! , 

?nmW,IJ " ^ UvWV (m - k)\k\{n - m + k)\ U n-m+k u k > 

* nm(q > X)= ^{wm) {n-k)\k\{ m -n + k)\ HkH --^ m>U - (8) 



Here flf are the Hermite polynomials of the following arguments, 
H+= Hi[i^h{n 2 q-ix)], Hf= H^^h^xq + ix)} 

[Hq{z) = l,H\(z) = 2z,Hi{z) = 4z 2 — 2,...], and C mn are the irrelevant for our purposes 
normalization factors. 

The wave functions © are normalized implying that the spectrum is pure pointd. However, 
for generic incommensurable frequencies, the spectrum is dense everywhere: however small the 
interval {E, E + dE} is, it contains an infinite number of eigenvalues, and this is true for any E. 
Obviously, this unusual property is related to the fact that the spectrum has no bottom. We 
want to emphasize, however, that, though this quantum system is unusual, it is not sick: the 
Hamiltonian © is Hermitian and the corresponding evolution operator is unitary. 

In the limit f^i — ► O2, the transformation ([3]) becomes singular. To understand what happens 
in this limit, it is best to look directly at the spectrum © and the wave functions (J7|), ©. 

If setting formally Oi = Q2 in equation (J0j(^ , we obtain 

E nm = D,(n- m) (9) 

meaning the infinite degeneracy of the spectrum at each level. The wave functions are reduced 
to 

" e -^V#+_ m5 m<n, 



(10) 

e qx H m _ n , m > n. 



We see that, at each level with a given n — m, an infinite dimensional Jordan block appears [5]. 



2 Mathematicians say that a self-adjoint operator H has pure point spectrum if all eigenfunctions of H belong 
to £2, forming a complete orthogonal basis there. 
3 See, however, the discussion below. 
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3 Jordan blocks and continuous spectrum. Unitarity 

The presence of Jordan blocks in the Hamiltoniar@ implies usually the loss of unitarity. Indeed, 
consider the matrix Hamiltonian 

*-(!!) <"» 

(it has only one eigenvector ip = ( ] with eigenvalue 1). It is straightforward to see that the 



x / 

time-dependent Schrodinger equation 

~ H * (12) 

has the following general solution, 

1 > —it , u I \ „-it 



m=a{ Q )* +b { £ )e-«. (13) 

When 6^0, the norm of ip(t) grows with time. The latter statement applies to the natural 
norm = iftip and also to any other positive definite norm ~ Mtp. If the norm is not 
positive definite, it is not true, as is clearly seen when choosing 



M 




1 



Such a norm is degenerate, however. It projects the full 2-dimensional Hilbert space where the 
Hamiltonian (jlip is defined to a one-dimensional subspace ip = f ^ ~\ . The dynamics in this 

subspace is unitary, the dynamics in full Hilbert space is not. 

This reasoning applies to any Jordan block of finite dimension and to any Hamiltonian inclu- 
ding such. However, in our case, the dimensions of the Jordan blocks are infinite. A remarkable 
fact is that in this case unitarity is not necessarily broken. In particular, a unitary evolution 
operator can well be defined for the free PU oscillator at equal frequencies. The novel feature 
compared to the inequal frequencies case is the appearance of continuous spectrum [31 [6] . 

Strong indications that the spectrum is continuous follow already from inspecting the spec- 
trum 0, the wave functions (|7|), and their fate in the limit Q\ — > fi2- 

• One can notice first of all that the discrete spectrum © is only obtained from © if setting 
f^i = while keeping n, m finite. It is possible to take the limit in such a way (this 
implies a kind of ultraviolet regularization), but a more natural approach is to allow n, m 
to be arbitrary large. When A = f2i — O2 is sent to zero and n, m are sent to infinity such 
that n — m is kept fixed and nA is kept finite, the energy ([6j) can acquire an arbitrary 
value and not only the discrete values Q . 

• In the limit A — > 0, the wave functions (J7]) are no longer normalizable, and this suggests 
that the spectrum is continuous. 

Let us now prove the continuity of the spectrum by constructing explicitly the wave functions 
of arbitrary energy E. At the first step, let us get rid of the terms tx x 2 and oc q 2 in the 



4 The points in the parameter space where such Jordan blocks emerge are called exceptional points, see e.g. [TT] 
and references therein. 
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Hamiltonian by representing an eigenfunction of ([2]) as ty(q,x) = e in xq (j)(q,x). The operator 
acting on cp(q, x) is 

H = \vl + xp q - n 2 qp x . (14) 

It is convenient then to perform the canonical transformation [3l [6] 

1 1 ] 



i 11 

Px -> Px, Pq -> Up q , X -> X + — p q , q -> —q + -T7y>Px- (15) 



giving the new Hamiltonian 

H' = ^^- + n(xp g -q Px ). (16) 

The transformation (|15p is the superposition of the scale transformational q — > g/O, p g — > f2p ? 
and the unitary transformation O — > e R Oe~ R with = ip x p q /(Afl). The eigenfunctions of the 
Hamiltonian (|14p is obtained from the eigenfunctions of the Hamiltonian (|16p aa^| 



<^,*)=exp{^^}^,.). 

The eigenfunctions of H' are known, 

<j>' lk (q',x;t) oc Ji(A;r) e -^e- 4 '( m+fc2 / 4 ), 

where (r, 0) are the polar coordinates in the plane (</, x) and / = 0, ±1, ±2, . . . are the eigenvalues 
of the angular momentum operator L = xp' q — q'p x - We obtain 



K exp {^^} [ J '(*V^ + nv) (^f) 

Introducing 

z = v / 0(x + iQq), w = iz = Vfi(ix + £lq), 
expanding the Bessel function and using the identity 



1/2 



e -it{m+k 2 /4) 



Adz 2 J 
we finally derive 

Hi,, nc r -it(m+k*/i) -itfgxY- (ik 2 \ m Hi +m (z)H m (w) 

m=0 x 7 v 

The structure of this expression is similar to that in equations (jlj), ([8]), but the meaning is 
different: the expansion goes now in the spectral parameter k 2 /O rather than in the parameter 
of the Hamiltonian A/yTJiTT^ as in equation 0. In addition, equation (|17p represents an 
infinite series rather than a finite sum. We can note that each level in the spectrum is infinitely 
degenerate. The eigenfuctions Vl/^, J k^+iVf e ^ c - nave the same energy. This is not a Jordan 



5 It amounts to a unitary transformation with S = e~^ qVq+VqqS> ln n,/2 , but this representation is not particularly 
useful. 

6 For a general theory of quantum canonical transformations see [12] and references therein. 
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degeneracy discussed above: the basis (|17p diagonalizes the Hamiltonian, and the functions 
with different I are distinguished by the eigenvalue of the "angular momentum" operator, 



L 



Xp q Q 1 

2n ~ 2 qPx + m 




3Q 3 q 2 



(we have rotated the operator L = xp' q — q'p x back to original variables), that commutes with 
the Hamiltonian ([2]). 

We have made two seemingly conflicting statements: (i) the spectrum of the Hamiltonian 
involves a set of infinite-dimensional Jordan blocks and (ii) the spectrum is continuous. How 
to reconcile them? To understand it better, consider the trivial free Hamiltonian 



H 



lcP_ 
2dx 2 



with the eigenfunctions 

, . r ik 2 

M/(x; t) = exp < ikx — t 



(18) 



Note now that not only (|18p . but also every term of its expansion in k, 



k° : 


^o(x 


t) 


= 1, 




k 1 : 




t) 


= x, 




k 2 : 




t) 


= t - 


ix 2 , 


k 3 : 


^ 3 (x 


t) 


= xt 


■ ^ 

IX 


k A : 


f 4 (x 


1) 


= t 2 


- 2itx 2 



~3~ 



(19) 



etc., satisfy the time-dependent Schrodinger equation (fT2l) . The functions (fT9l) represent polyno- 
mials in x and, what is especially noteworthy, in t, much similar to non-stationary solutions (|13p 
characteristic to Jordan block Hamiltonians. 

On the other hand, the functions (|19p grow at large x, are not normalizable, and do not form 
a basis of a reasonable Hilbert space. The only way to define the latter is to sum over all ^ n {x\ t) 
with a proper weight and go back to the standard continuous spectrum wave functions (I18p . 
which can be dealt with by introducing a box of finite length L (where the definition of Hilbert 
space presents no problem) and sending then L to infinity. 

By the same token, individual terms of the expansion of (|17p in k represent nonstationary 
solutions of the Schrodinger equation with the Hamiltonian ([2]). For 1 = 0, these solutions are 



. e -ifl 2 qx 



[t - i[x" + \L"q") \t 
t 2 - 2it(x 2 + n 2 q 



k 

k A : 



2 + n 2 q 2 )]e- in '^ 

2 . (n 2 q 2 +x 2 ) 2 



iqx + 



sn 2 



-ifl qx 



(20) 



etc. The first function in (I20p coincides with the wave function (I10p with n — m = 0. The higher 
functions may be interpreted as its non-stationary "descendants" associated with the presence 
of the infinite-dimensional Jordan block at the level E nm = 0. 

However, as was also the case in the previous example, these descendants grow at large x, q 
and only their superpositions (fTT]) form a proper Hilbert space basis. 
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4 Including interactions 

We have learned in the previous section that, in spite of unusual features (the absence of the 
ground state), the free PU oscillator represents a benign unitary quantum system both when 
fii 7^ O2 and when Oi = f^- It is especially clear in the inequal frequencies case when the 
Hamiltonian can be brought in the form ([5]). When two different oscillators are decoupled, it 
does not matter much that the energies of one of the oscillators are counted with the negative 
sign. Each oscillator lives its own life and the energy sign is basically a bookkeeping issue. 

The danger arises when the oscillators start to interact. The subsystem 1 can give energy 
to subsystem 2 such that the absolute values of the energies of both subsystems increase. This 
brings about a potential instability that may lead to collapse and loss of unitarity. We will 
see that, generically, such collapse occurs, indeed. But not always. When the interaction 
Hamiltonian has some special form, there is no collapse and unitarity is preserved. 

Consider the Hamiltonian (acting on 4>(q,x)) 

H = ~\^ + ~ ix J- + «« 4 + fa 2 * 2 + 7* 4 - (21) 

Let us find out first whether the spectrum is still unbounded from below as for the free PU 
oscillator. The answer is — yes, it is. To see that, choose the variational Ansatz 



f Aq 2 „ Cx 2 ) [AC , . 

x)=exp<^ iBxq ' ' 

The variational energy is 

F(A B D - C + B2 + BQ2 * 3a 37 13 



4 AA 2A 2C AA 2 AC 2 AAC 

By varying A, B, C, one can make it arbitrarily negative. Indeed, choose first A very large such 
that 

Keeping C finite and increasing B, we can make —E arbitrary large. In [5], we used a different 
variational Ansatz with which the energy of the ground state could not be made arbitrarily 
negative. However, as the energy can be lowered indefinitely with a more general Ansatz (|22p . 
we conclude that the ground state is absent. 

When a or (3 or 7 are nonzero, the classical dynamics of the Hamiltonian (|2ip involves collapse. 
In the certain region of parameters, there is an island of stability around the perturbative 
vacuum [2], defined as the point q = q = q = q^ = 0. With initial conditions at the vicinity 
of this point, the trajectories display benign behaviour peacefully oscillating near the origin. 
Generically, however, the trajectories go astray reaching infinity in finite time. 

The simplest well-known system displaying such behaviour describes 3D motion of the particle 
with the attractive potential V(r) = —n/r 2 . Classically, for certain initial conditions, the particle 
falls to the center in finite time. The quantum dynamics of the system depends in this case on the 
value of k. If mK < 1/8, the ground state exists and unitarity is preserved. On the other hand, 
if mK > 1/8, the spectrum is not bounded from below and, which is worse, the spectral problem 
cannot be well posed until some extra boundary conditions are set at the vicinity of the origin. 
The spectrum depends then on these boundary conditions [13] . Setting such boundary conditions 
is tantamount to regularizing the singularity of the potential (the spectrum thus depends on the 
way it is done). Without such regularization, the probability leaks out through the singularity 
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Figure 1. A typical classical trajectory of the Hamiltonian (|24|) . 



and unitarity is violated. One can conjecture that the presence of classical collapsing trajectories 
and the absence of the ground state always imply violation of unitarity. We do not know whether 
a mathematical theorem to this effect can be proven, but, heuristically, it is suggestive. Indeed, 
if the spectrum involves the states with arbitrary low energies, the corresponding wave functions 
should have main support near the singularity, where theory needs regularization to prevent the 
leakage of probability. 

In other words, the nonlinear PU oscillator is in most cases sick: collapsing classical trajec- 
tories lead to the quantum collapse. However, in some cases there is no collapse. Consider the 
PU oscillator with inequal frequencies and let us add to the Hamiltonian ([5]) the potential 



V(X u X 2 ) = X(X 1 -X 2 )(X 1 +X 2 f, A>0. 



(23) 



The numerical solution to the classical equations of motion, 
dV 



x 1 + n(x x + 

x 2 + Q 2 X 2 



dX l 
dV 
~dX 2 



0, 
0, 



do not display any collapse. The amplitude of oscillations somewhat grows with time, but this 
growth is rather smooth (see Fig. [I]). 
The Hamiltonian 



H 



+ \{X 1 -X 2 ){X l +X 2 ) 



(24) 



is a close relative to the super symmetric Hamiltonian considered in The bosonic part of the 
latter is 



H = P P + D{n 2 x + Xx 3 



(25) 



where (p,P) are the canonical momenta of (x,D). By introducing X\ )2 = {x ± D)/y2, it is 
reduced to the Hamiltonian (|24p with f2i = = The system ()25[) is integrable and the 



7 To avoid confusion, it is worth reminding that the limit Oi — > is singular and the Hamiltonian ([5]) with 
Qi = f^2 does not describe the PU oscillator with equal frequencies. 
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E<0 



EXI 



Figure 2. Stokes wedges. 

classical trajectories are expressed analytically into Jacobi elliptic functions. In this case, the 
amplitude of oscillations grows linearly with time, but this does not represent a disaster as 
the trajectories do not run into infinity in a finite time. The Schrodinger equation for the 
Hamiltonian (|25h can be solved analytically [4]. The spectrum is continuous, with eigenvalues 
lying in two intervals ] — oo, — 0] U [Q,oo[ plus the eigenvalue {0}. The quantum spectrum 
for the Hamiltonian (|24p should have the same qualitative features. Besides (|23p . there are 
other interaction potentials providing for the same behaviour. One of them is V(X\, X2) = 
X(Xi — X2) 3 (Xi +X2). Note that, being expressed in terms of the original variables x, q, p x , p q , 
the Hamiltonian (|24|) acquires a rather complicated form and the corresponding Lagrangian 
depending on q and its time derivatives is even more complicated. Natural nonlinear extensions 
of the Lagrangian ([T]) seem all to involve classical and quantum collapse. 

5 Complexification 

The spectral Sturm-Liouville problem is defined when the operator in question (Hamiltonian) 
is given and boundary conditions are specified. It is possible thus to have different spectral 
problems associated with a given Hamiltonian. 
Take a simple oscillator, 

h = ^(p 2 + n 2 x 2 ). 

If considering the spectral problem at real X and requiring that wave functions vanish at infinity, 
we obtain the usual spectrum 

E n = fi(n + l/2). 

The wave functions, \P cx e~ x2 ^ 2 , continued analytically to complex X, vanish asymptotically 
not only on the real axis, but in the whole Stokes wedge with |arg(X)| < 7r/4. Nothing prevents, 
however, to consider another spectral problem and require that wave functions vanish at large 
imaginary X. The corresponding eigenfunctions are ^ n cx H n (iX)e x I 2 . They vanish in the 
complementary Stokes wedge as shown in Fig. [2j The spectrum is now 

E n = -«(n+l/2). 

We have seen that the free PU oscillator with different frequencies is equivalent to the combi- 
nation of two oscillators ([5]) with the spectrum ([6]). But the latter is true only if q, x, p q , p x and 
hence X\p are assumed to be real. One can on the other hand assume that, while X\ is kept 
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real, X 2 is imaginary. In this case, the sign of the energy eigenvalues of the second oscillator is 
reversed, and one obtains a positive definite spectrum 

E nm = ( n + - ] fii + ( m + - ) n 2 , n,m = 0, 1, ... (26) 



with no trace of ghosts. Another way to obtain this result is the following El [9]. Let us consider 
the Hamiltonian ([2]) and assume x being real and q being imaginary therqj, q = iy, p q = —ip y - 
The Hamiltonian acquires the form 

H = Y~ lxp y + 2 + 2 (27) 

The second term in this expression is complex which may lead to worries whether the Hamilto- 
nian (|27p makes sense. This worries are not grounded, however. Though the Hamiltonian (|27p 
is not manifestly Hermitian, it belongs to the class of crypto- Hermitian Hamiltonians with real 
spectrum |14| . The Hamiltonian (|27p can be diagonalized to the manifestly Hermitian form, 



-Hdiag - ^ 1 ^ • l^oj 

by a complex canonical transformation (a brother of the transformation ([l])), 
_ 1 Ofx - zp y Slfy - ip x 



P — ~ Q P' x ^ p — p y + i&2 x ^29) 

This transformation amounts to a certain similarity nonunitary transformation, 

#dia g = e~ s e- R He R e s (30) 



with 



■ PxPy xyQiQ 2 \ ln + ^2 



20i^ 2 2 J Oi-0 2 

and 

5 = zll ^ ni (yp y +p y y)- 

To find the eigenfunctions of (|27|) . one can either rotate away the eigenfunctions of 
^ = e R e 5 ^diag or proceed in the same way as in [5] representing the wave functions as 

*nm = <?Wexp j- — Li_^(a; 2 + Oif^y 2 ) - ftifi 2 a;y 
The operator acting on (j) nm (x,y) is 

6 = "55? - 4y + W"' + ">> + + ^TT- 



8 From the viewpoint of the original Lagrangian fl}, it is not so natural to assume that the coordinate q is 
imaginary, while its time derivative x — q is real. But, formally, one is allowed to do it. 
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Introducing 

z = \f^h{x + Q 2 y), 
it acquires the form 



w 



= fii 



1 

Hlh 2 



d_ 

' dz 



+ n 2 



tt 2 (x + Qiy), 



1 d 2 d_ 

2 dw 2 dw 



(31) 



The eigenfunctions of O are polynomials. When the dependence on the variable w or on the 
variable z is suppressed, they are conventional Hermite polynomials, 



>n0 



H n (z) 



H m (w). 



When both n ^ and m ^ 0, the eigenfunction is 



4>nm(x,y) 



min(n,m) 

£ ' 

k=0 



m\(n — to)! 



»1 + ^2 

4 v / n^y (to — k)\k\{n — to + A;)! 



—m+k 



(32) 



The expression (|32p is very similar to ([5]) and is obtained from the latter by substituting iA = 
— — ► — + an d replacing the arguments i\/TF[(Q2Q — ix) — > z, ^y^l2(0,lq + ix) — > w. 
There is a price, however, that one has to pay for getting rid of negative energies in the 
spectrum. The Hamiltonian f|2Tj) is not Hermitian and this means that the conventionally 
defined norm = JJ dxdy \*B(x, y)\ is not preserved during evolution. In addition, the 



eigenfunctions (|32p do not represent an orthogonal basis with respect to this norm. True, 
crypto-Hermiticity of the Hamiltonian implies that a norm with respect to which the evolution 
is unitary can be defined, but this norm, 



1^1 



(q* e - 2R y 



has a complicated nonlocal structure. 

The transformations (|29p are singular at f^i = f^- This suggests that, as it was the case for 
the standard PU oscillator, the point f^i = is exceptional also in the nonstandard realization. 
The emergence of Jordan blocks in this limit was noticed back in [7], but the best way to see 
that is to follow, as Bender and Mannheim did [9], the approach of [5] and to explore the fate 
of eigenfunctions (|32|) in the limit f^i — > 

In this limit, the variables z and w coincide. The eigenfunctions (I32h for a given n + m also 
all coincide in this limit, (j) nm (z) ~ H n+m (z). To derive the latter, note that the operator (|3ip 
acquires in the limit J7i — > £l 2 the form 



On 1= Q 2 = ^ 



a 



d 2 n o ; 

+ 2z— + 1 



dz 



and its eigenfunctions are simple Hermite polynomials, indeed. After massaging equation ([32 
a little bit, we derive, as a byproduct, a nice mathematical identity 



min(n,m) 

E (- 2 ) j 

3=0 



n\m\ 



j\(n - j)\{m - j)\ 



H n _j(z)H m _j(z). 



In contrast to the standard PU oscillator where the Jordan blocks had infinite dimension, their 
dimension is finite here. We have a single vacuum state with energy E = the Jordan block of 
dimension 2 at the level E = 2Q, the Jordan block of dimension 3 at the level E = 30, etc. At 
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each level, there is a finite number of different nonstationary solutions to the time-dependent 
Schrodinger equation and one cannot construct, as we did before, bounded in x, t combinations 
that have the meaning of continuous spectrum wave functions. As a result, unitarity is violated^. 

As was mentioned, the real problem associated with the ghosts is the quantum and/or clas- 
sical collapse, which might appear only in interacting theory. It is not so easy to include 
interactions and analyze their effects in the Bender-Mannheim approach. First, it is not clear 
whether the Hamiltonian (|27p with the interaction term like Ay 4 is still crypto-Hermitian. The 
canonical transformations (|29p [or, which is equivalent, the similarity transformation (I30j) ] kill 
the complexity —ixp y , but bring about new complexities coming from the interaction term oc y 4 . 
These complexities have the form AH^g oc iA(O^PiA| — PfX2). We failed to find a modified 
canonical transformation that would kill all complexities. 

Even if one succeeds in finding a nonlinear crypto-Hermitian generalization of the Hamilto- 
nian (I27p . it would be difficult to analyze its spectrum and to find out whether the ground state 
is still present there. In particular, it would be difficult to implement variational estimates due 
to a complicated nonlocal norm. 



6 Conclusions 

Our main message is that the ghosts (the absence of the ground state in the spectrum), do 
not represent a serious problem for free theory. Thus, it is not necessary to cope with them 
there. The spectrum of the free PU oscillator runs from -co to +00 representing a pure point 
spectrum that is dense everywhere when 7^ Q2 and a continuous spectrum when = ^2- 
The latter can also be interpreted via infinite-dimensional Jordan blocks that appear in the limit 
^1 —> ^2 [5]- In spite of the absence of the ground state, unitarity is preserved. 

Generically, interactions break unitarity due to quantum collapse phenomenon. However, 
there are certain special cases when the theory involves neither classical nor quantum collapse. 
A particularly interesting example is the exactly soluble nonlinear Hamiltonian (125j) that appears 
in the context of super symmetric higher derivative quantum mechanics [^j^l . 

By analytical continuation to complex coordinates, it is possible to consider a nonstandard 
realization of the free PU Hamiltonian with the positive definite spectrum (|26p . When f2i ^ Q2, 
a unitary evolution operator can be defined. The presence of the ground state may be considered 
as a certain advantage of this realization compared to the standard one. However, 

1. It is misleading to say that choosing the nonstandard realization "solves the ghost prob- 
lem". The analysis has been performed only for free theory and, in the free case, there is 
no serious problem to solve anyway. 

2. In contrast to the standard one, the nonstandard PU oscillator at equal frequencies is not 
unitary. 

3. As far as the original problem with the Lagrangian (pQ) is concerned, the nonstandard 
realization does not look natural: (i) it is not so natural to assume that the coordinate q 

9 In [9], a certain norm was constructed that is conserved during evolution. As a result, Bender and Mannheim 
claimed that the system with the Hamiltonian (|27[) in the equal frequency limit is unitary. However, their norm is 
nilpotent and, as was discussed at the beginning of Section 02 it amounts to projecting the original Hilbert space 
onto a complicated unnaturally defined subspace. Even though the dynamics in the latter is unitary (and the 
Hamiltonian loses its Jordan block structure and is Hermitian), this is not the unitarity and Hermiticity in the 
standard meaning of these words. 

10 The Hamiltonian (|25p had continuous spectrum, like the free PU oscillator with equal frequencies. Another 
interesting nonlinear higher-derivative Hamiltonian, 

H = pP + D(Q. 2 x + \x 3 ) - |(D 2 + P 2 ), 

has the pure point spectrum that is dense everywhere, like the PU oscillator at different frequencies. 
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is imaginary while its time derivative x = q is real; (ii) the norm in the space (y = — iq, x) 
that is preserved during evolution has a complicated nonlocal structure. 

4. It is not known yet what happens in the framework of nonstandard realization when 
interactions are included. 
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